CASE FILE 
COPY 

NASA TECHNICAL 


CNJ 

CS| 



on 


ANALYSIS OF HYPERYELOCITY 
PERFORATION OF A VISCO -PLASTIC 
SOLID INCLUDING THE EFFECTS OF 
TARGET-MATERIAL YIELD STRENGTH 


by Robert G. Thomson 

Langley Research Center 
Langley Station, Hampton, Va. 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • APRIL 1965 




NASA TR R-221 


ANALYSIS OF HYPERVELOCITY PERFORATION 

OF A VISCO-PLASTIC SOLID INCLUDING THE EFFECTS 

OF TARGET- MATERIAL YIELD STRENGTH 

By Robert G. Thomson 

Langley Research Center 
Langley Station, Hampton, Va. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For solo by the Office of Technical Services, Deportment of Commerce, 
Wash i ngton, D. C . 20230 -- Price $2.00 




ANALYSIS OF HYPERVELOCITY PERFORATION 


OF A VISCO-PLASTIC SOLID INCLUDING THE EFFECTS 

OF TARGET -MATERIAL YIELD STRENGTH 

By Robert G. Thomson 
Langley Research Center 


SUMMARY 


A visco-plastic flow theory for the solution of hypervelocity perforation 
of thin plates has been investigated. A short-time analytical solution of the 
governing equation including the effects of yield strength of the target mate- 
rial was obtained and expressions were derived for the velocities, displace- 
ments, strains, and strain rates present in the target material during hyper- 
velocity impact. In addition, pertinent parameters have been varied and a 
comparison has been made with the solution in which the yield strength was neg- 
lected. The results indicate that the velocities, displacements, and stresses 
determined from the present solution are at times markedly different from the 
results obtained when the yield strength is neglected. The effects on the radius 
of perforation of the plate material when the yield strength is included were 
found to be not as significant as those of the solution neglecting target yield 
strength and resulted in little variation in the computed perforation radii. 


INTRODUCTION 


A visco-plastic flow theory was first proposed for the hypervelocity per- 
foration of thin plates by rigid cylindrical projectiles by F. A. Bakhshiyan in 
194-8. (See ref. 1.) In this analysis a circular cylindrical projectile was 
considered to impact upon a thin target plate and the resulting radial visco- 
plastic boundary was treated as time dependent. The solution presented, however, 
was limited to the special case of an infinite mass projectile. In references 2 
and 3 an analytical solution was obtained for a problem similar to that of ref- 
erence 1 except that the yield strength of the target plate was taken to be zero 
(that is, the visco-plastic boundary was taken to be at infinity) and the pro- 
jectile mass was treated as finite. The yield strength of the target material 
was introduced later in this analysis to establish a separation criterion for 
perforation. An approximate two-dimensional analysis of the impact problem as 
stated in reference 2 is presented in reference 4*. This solution indicates 
that the radial velocities are about one order of magnitude less than the axial 
velocities and thus helps to justify the use of a one -dimensional approach. 



An extension of the work presented in references 2 and 3 is presented in 
reference 5 in which an attempt is made to include the effects of target mate- 
rial yield strength in the governing differential equation and associated 
boundary conditions. The results of the analysis are given in finite series 
form and a parametric study is presented in which the effects of variations in 
the pertinent parameters on the radius of perforation are determined. The 
solution as given in reference 5, however, does indicate that for the case of 
the infinite mass projectile, the velocity of the projectile decreases slightly 
after impact. This result is probably due to the fact that only a finite num- 
ber of terms from an infinite series was used. 

In the present analysis an analytical solution to the problem of a rigid 
projectile of finite mass impacting upon an infinite plate is obtained in which 
yield-strength considerations have been included in both the analysis and in 
the determination of the separation criterion. The solution obtained in the 
present report differs from that reported in reference 5 even though the basic 
problem is identical. The present solution agrees in the limiting case when 
the yield strength is set equal to zero with the solution presented in refer- 
ences 2 and 3 in which yield strength was neglected. In addition, the expres- 
sion for velocity will limit, in the case of the infinite mass projectile, to 
the condition that the velocity of the projectile remains constant after impact. 

The present report also shows how the effects of including yield strength 
in the analysis compare with the results obtained for the same conditions when 
yield strength is neglected. Some pertinent parameters such as initial velocity 
and mass of the impacting projectile as well as the plate thickness and viscos- 
ity are varied in order to study their effects on the displacement, velocity, 
strain, and strain-rate distributions, and on the radius of perforation as com- 
puted with and without target yield strength. 


SYMBOLS 


A,B coefficients of Bessel functions 

a radius of projectile 

Cp dynamic ultimate yield strain in shear 

erf error function 

erf' derivative of error function 

erfc complementary error function 

g Q initial projectile velocity 


2 



h 


plate thickness 


5 

H 


2:rtah 

M 

2jta^hp B 

M ’ 


also used for Hankel function 


i n erfc x 
In 


k 

K 

K n 

Mp 

M 


nth integral of erfc x (see appendix, eq. (A 38 )) 
modified Bessel function of first kind of order n 


Bessel function of first kind of order n 


dynamic target yield stress in shear 
ka 


Bingham number. 


Vo^ 


modified Bessel function of second kind of order n 
mass of projectile 

mass of projectile and plug of plate material of radius 

Mp + jta^hp 


n,p integers 

r radial distance 

r nondimensional radial distance, — 

s Laplace transform parameter 

Sp Struve function of order p 

t time 

t nondimensional time, —t 

a^ 

V axial velocity 

V transform of V 

9 angle in plane of target (see fig. l) 



Vo 

v_ 

Vo 

dV _a_ 

Sr V D 


initial velocity of projectile after impact 
nondimens ional velocity 

nondimens ional shear strain rate 
axial displacement 



dw v 
Sr aV 0 

Yn 


7 

v 


nondimensional displacement 

nondimensional shear strain 

Bessel function of second kind of order n 

shear strain, ^ 
dr 

dynamic coefficient of viscosity 
kinematic coefficient of viscosity, ^ 


p mass density of plate material 

r rz transverse shear stress 

z axial coordinate 

Subscripts: 
cr critical 

p perforation 

A dot over a symbol denotes differentiation with respect to time t. 


ANALYSIS 


Governing Equations 

In the present analysis a rigid cylindrical projectile is considered to 
impact upon a thin infinite plate. The resulting perforation of the plate by 
the projectile is assumed to be a simple shear-plugging perforation in which 
only the transverse shear stresses act to resist the inertia of the impacting 
projectile. The perforation is also considered to be axially symmetric and the 
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shear stress is taken to he constant through 
the thickness of the plate. The resulting 
deformations w of the plate are then 
represented as functions only of the radial 
coordinate r and the time t, both being 
independent of the axial coordinate z and 
the circumferential coordinate 0. (See 
fig. 1(a).) Hence, by taking the sum of the 
forces in the axial direction on a circular 
ring element of plate material ( see 
fig. 1(b)), the basic equation of motion 
for this simple shear perforation model can 
be written as 



(a) Finite rigid projectile 
and infinite plate. 


^T- 


rz 


Sr 



( 1 ) 



( T r z d T r z dr) h^r + dr)d 


The plate material considered herein is 
assumed to behave like an incompressible, 
visco-plastic Bingham solid in which the 
deformation accompanying transverse shearing 
commences only after the yield strength of 
the material has been reached. (See 
ref. 6.) When the value of the transverse 
shear stress falls below the yield stress 
of the material (or when the rate of defor- 
mation becomes equal to zero), visco- 
plastic flow ceases and the material is 

assumed to be rigid. The relation between the shear strain rate and shearing 
stress for the case of simple shear perforation can be written as (see ref. 6) 


(b) Circular ring element with 
equilibrium forces. 

Figure 1.- Geometry and coordi- 
nate system of perforation. 


T rz 




(2a) 


or, since the 


sign 


dw 

dr 


is always negative in this case, 


r 


rz 



(2b) 


where 

|i dynamic viscosity of target material 

k dynamic yield stress in shear of target material 

= shear strain rate 

or or 
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w, and strain y, 


and the relations between the velocity V, axial displacement 
are 

v = dw. 

St 



by _ SV 
St Sr 


Substitution of equation (2b) into equation (l) yields the governing linear 
differential equation 




S 2 V 1 Sv 1 SV _ 1 k 

^2 + rSr~vSt“rn 



(3) 

where v 
are taken 

is the 
to be 

coefficient of kinematic viscosity. 

The initial conditions 




V = 0 

(t 

=0 ; r > a) 

(*0 



V = V 
v v o 

(t 

= 0, r S a) 

(5) 

where V 0 

is the 

initial velocity of the projectile and plug 

combination. 

In 


this study V Q was determined from the conservation of momentum at the instant 
of impact in which the projectile-plug combination is assumed to be rigid and 
to act as a unit. Thus 


V 


o &Q 


% 

M 


(6) 


where 

g Q initial velocity of free projectile 

Mi mass of projectile 

M = Mi + jta 2 hp 

The boundary conditions are 


-Ck + ^ 


S V 
Sr 



(t > 0, r = a) (7) 
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and 


V -» 0 


(t ^ 0^ r -> oo ) (8) 


where 


t - 2rtah 

^ M 


Method of Solution and Resulting Equations 

An analytical solution is obtained in the present paper by Laplace trans- 
form techniques similar to the approach used in reference 5* The homogeneous 
solution can be written directly in terms of modified Bessel functions. The 
particular solution, however, is obtained by utilizing Struve functions of order 
zero rather than by the method of variation of parameters used in reference 5* 
The solution to the governing differential equation (eq. (3)) and its associated 
boundary conditions (eqs. (4), (5)* (7)j> and (8)) is discussed in detail in the 
appendix. The solutions obtained are "short time" solutions in t corre- 

2 

sponding to values of the transform parameter s > — or t < (see 

appendix) resulting from the use of asymptotic expansions in the transformed 
state. It is shown in the appendix that the velocity (eq. (A3^))» displacement 
(eq. (A37) ), shear strain (eq. (A 36 )), and shear strain rate (eq. (A35)) can be 
expressed in dimensionless form in terms of a dimensionless radius r and a 
dimensionless time t as 


_v_ 

v 


l 

ft 



r ~ 1 + 4 k( 1 - H)ti 2 erfc + l6KHt 2 i^erfc E- S -l 

2 ft 2^1 2ft 


+ Rn 


2ft i erf c A + 8K(l - H) (t)^ 2 i 5 erfc r . ~ 1 + 32KH(t)^ 2 i 5 erf c r ~ 1 

2ft 2/l 2/l 


R 2 


4ti 2 erfc r ~ 1 + l6K(l - H)t 2 i^erfc r ~ - 1 + 64KHt 5 i 6 erfc r 

2 Jf 2{? 2ff 


((r - X) £ 0) (9) 


^ (i 2 erfc r ~ 1 + 4X(l - H)ti i *‘erfc r - ~ 1 + lfiKHt^erfc I— i 
a V o ft I 2ff 2^T 


2\[t i^erfc 



8K(l - H)(t) 5 / 2 i 5 erfc r ~ 1 + }2KH( t ) 5 / 2 i T erfc 

2ff 



+ R. 


4 

4ti erfc 


1 I + 16k( 1 - H)t 2 i 6 erfc r " 1 + 64KHt 3 i 8 erfc r ~ 1 . 

2/l 2^f 2 /? 



2r 


((r - 1) * o) (10) 
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Sv V 

Sr aV Q 


SV _a_ 
Sr V 0 


where 


-p<- erf c - -~ - J; + R^i 2 erfc - — — i + 2 /tll?^ - K(l - H)| l^erfc 1,-1 

yr I 2|ft 2|t 2^1 L “ J 2|/r 


l+t[i 5 + K(l - H)R 3 ] i^erfc + 8 (t) 3 / 2 [k(1 - H)R 4 - Kh] i 5 erfc 

l j ojt l j 


r - 1 




+ 16 ? 


|kHR 5 + K(1 - H)R^ll 6 erfc — '? ■ + 32(t) 5 / 2 KHR4i T erfc r ~ 1 

L ” J 2yt 2\fi 


+ 6^ %HRci®erfc E . - 3a . 

5 2 Jf 


k? 

2 ? 


((r - 1) > 0 ) (11) 


(F-l ) 2 

1 A 1. 


+ PU erfc — ~ ^ + 




2||t £4 - K(l - H)1 i erfc 

2\ft 


+ Vtfe + K(l - H)R^j i 2 erfc r ~ 1 + 8 (t) 5 / 2 |K(l - H)R 4 - Knji^erfc 
+ l 6 ? fjCHRj + K(l - H)R^1 i^erfc r ' — — + 32 (t) 5 / 2 jQjR^ 5 er f c 

L J 2 \jt 


r - 1 

2^jt 


r - 1 






r - 1 

■> tE i 

((r - 1) > 0 ) 

2 ft 


V / 


H = 


2 rta^hp 


IT 


M = Jta^hp + M-^ 
ka 


K = 


vVq 


t =^t 

a^ 


r = ^ 
a 
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r 2 = 


2-r - - A + H[i - 2 1 + H2 


128f 2 64F 128 \8r 8 


R* = - 


5_ _ 1 

8r 8 


+ H 


Ri, 


= + -2- + 2(2 + 5^ _ 

128F 2 6^r 128 8lr ^ 1 


Rr = 


. * t 5 . + _1_ + -1_ + JLfe . JJ 

256?5 128 r 2 256 r l6\r r 2 


Hf. 

2r 


and i n erfc r ~ r is the nth integral of erfc r ~ (See eq. (A58) • ) 


!/? 




It can be seen from these equations that if the yield strength k is 
neglected, that is, if K = 0, the expressions for velocity, displacement, 

shear strain, and shear strain rate reduce to expressions identical to those 
presented in reference 3* Also in the case of the infinite mass projectile in 
which H = 0 (for = oo^, the expression given in equation (9) fen velocity 

is shown to reduce to V = V Q at r = 1 for all values of t . 

Typical examples of the radial distributions of strain and strain rate as 

expressed in equations (ll) and (12) are shown in figures 2 and 3* Typical 

velocity and displacement distributions as given by equations ( 9 ) and (lO) are 
plotted in subsequent figures (figs. 9 and 10). 


Separation Criteria 

In order to determine the radius of the plug of plate material which is 
sheared (or perforated) from the plate upon impact, separation criteria were 
established. The criteria are based on the assumption that separation of the 
projectile and plug of plate material from the plate occurs when the plate 
material can no longer transmit shear stress. If the plate material is con- 
sidered viscous, it can transmit high shear stress even though visco-plastic 
flow is occurring. This condition of high shear stress and visco-plastic flow 
exists immediately after impact when the strain rate is at a maximum and the 

rW 

viscous stress p—L is much greater than the yield stress k of the material, 
dr 

(See eq. (2b).) The plate material is considered to be highly viscous when the 
strain rate is greater than k/p but is considered susceptible to separation 
once the strain rate falls below k/p. In other words, the strain-rate 
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Figure 2 .- Strain as a function of radius for h = 0.25 inch ( 0.635 cm); 

Ml = 3.56 X 10-6 slug ( 51.95 mg); (i = 100 - SlUg (o.klffi £§_\j 

ft-sec ^ cm?/ 

■ 2 Sf l 2 - 68 A) ; k = 100 ' 000 psi 

68.95 J2LV a = 0.01+688 inch (0.119 cm). 

, cm^y 


So = 3,940 fps (l ,201 a); p = 5 


criterion for separation is assumed to be 


(it ■ f (15) 

A second criterion deals with the shear strain of the material which at the 
moment of impact is zero and increases thereafter. In order for the plug to 
separate from the plate, the material not only has to be considered as having a 

sufficiently small strain rate but also must contain large shear 

strains as well. The material is assumed to contain large shear strains when 
the shear strain exceeds a certain magnitude. The magnitude chosen in this case 
is the dynamic value of the ultimate shear strain of the material. Thus, the 
second criterion can be written as 


■ (it ■ (i1 

where Cp is the dynamic ultimate strain in shear. (A value of Cp of 0.02 
was assumed in this paper for both aluminum and steel.) For a further discus- 
sion of these separation conditions, see references 2 and 3 . 
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Figure 3*- Strain rate as a function of radius for h =0.25 inch 

(0.635 cm); M X = 3.56 X 1CT 6 slug (51-95 mg); n » 100 Ug - 

ft- sec 

(0.14-786 go = 3,9^0 fps (l,201 |); P = 5-2 (2.68 

k = 100,000 psi (68.95 — a = 0.04688 inch (0.119 cm). 

\ cmv 

The radius of separation of the plug material and the time of separation 
can most readily be determined by a graphical procedure. In figure 2 is shown 
a typical nondimens ional set of curves of strain as a function of radius and 
in figure 3 of strain rate as a function of radius for a projectile of mass 
3.56 X 10 slug (51*95 mg) and initial velocity 3*9^0 fps (l,201 m/s) impacting 
on a l/l+ inch (0.635 cm) thick aluminum plate. (Conversion factors for the 
units used herein are given in table I.) A plot of the time parameter t as a 



function of r can nov be obtained for critical strain from figure 2 by con- 
structing a horizontal line representing the critical strain (eq. (l4)) and 
determining the intersections with the strain curves. In a similar manner the 
critical strain-rate curve of t as a function of r can be determined from 
figure 5 . The radius of perforation and the time of separation are then deter- 
mined from the intersection of the critical strain and critical strain-rate 
curves as plotted in figure 4. Thus, for this particular case the radius of 
perforation and the time of separation are rp = 1.740 and t = 0 . 250 , 
respectively. 


TABLE I.- CONVERSION FACTORS FOR UNITS 


Parameter 

U.S. Customary Unit 

SI system 

Conversion factor 
(*) 

Mass 

slug 

milligram 

lb . 6 x 10 6 

Velocity 

fps 

meters 

second 

0. 5048 

Length 

inch 

centimeter 

2.54 

Coefficient of viscosity 

lb-sec 

ft 2 

newton- sec 

centimeter^ 

0.004788 

Yield strength 

psi 

kilonevton 

centimeter^ 

6.895 x 10‘ 4 

Density 

slug 

f t 3 

grams 

centimeter^ 

0.51558 


"’‘Multiply value given in U.S. customary unit by conversion factor to obtain equiv- 
alent value in SI unit. 



Figure 4.- Graphical representation for determining perforation 
radius and separation time, h =0.25 inch (0.635 cm); 

Ml = 3.56 X 10-6 slug (51.95 mg); p = 100 SlUg /0.4788 — 

ft- sec \ cm 2 / 

go = 3,940 fps (1,201 p = 5.2 (2.68 JLj; 

k = 100,000 psl (68.95 a = 0.04688 Inch (0.119 cm). 

\ cm 2 / 
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Note that four regions are indicated in figure 4. In these regions the 
predominant material behavior could be described as containing: large strain 

rates and large strain (region I); large strain rates and small strain 
(region II); small strain rates and small strain (region III); and small strain 
rates and large strain (region IV). Only in region IV are both critical compo- 
nents of the separation criteria satisfied; therefore, it is assumed that only 
in this region does separation occur. 


EFFECTS OF TARGET -MATERIAL YIELD STRENGTH 
ON PERFORATION RADIUS 


Expressions are derived in the appendix and presented in the previous sec- 
tion for the velocity, displacement, shear strain, and shear strain rate written 
in nondimens ional form* These expressions contain pertinent parameters associ- 
ated with the perforation problem such as projectile mass, initial velocity, 
target density, viscosity, yield strength, and thickness. Variations in these 
pertinent parameters are now investigated to determine their effects on the 
perforation radius. Results are obtained, with the aid of electronic computing 
machines, for the present solution which includes the effects of yield strength 
and for the simplified solution in which yield strength is neglected (that is, 
in which K is set equal to zero in eqs. (9)* (lO), (ll)j and (12)). 


Perforation Radius as a Function of Initial Projectile Velocity 


Aluminum plate .- In figure 5 is shown the variation in perforation radius 
with initial projectile velocity for a l/4 inch ( 0.635 cm) thick aluminum plate 
being impacted by a 3/64 inch ( 0.119 cm) radius rigid cylindrical projectile. 

The mass of the projectile was taken to be 3-56 X 10 "6 slug (51*95 mg) .which for 
aluminum would correspond to a projectile 0.0642 inch (O.I 63 cm) long. The 

dynamic coefficient of viscosity of the plate is 100 ^" s ^ c ( 0.4788 and the 

/ IcN \ ft 2 \ cm 

dynamic yield strength 100,000 psi (68.95 V As can be seen from figure 5> 

\ cm 2 / 

including material yield strength in the analysis has little effect on the per- 
foration radius Tp for the range of velocities shown. (The maximum difference 
between the case K / 0 and the case K = 0 is of the order of 5 percent.) A 
calculation was made for an initial projectile velocity of 778 fps (237 m/s) and 
the projectile did not perforate the plate. 


Steel plate .- In figure 6 is shown the variation in perforation radius with 
initial projectile velocity for a 0.088 inch ( 0.2235 cm) thick steel plate being 
impacted by a 3/64 inch (0.119 cm) radius rigid cylindrical projectile. The 
mass of the projectile was taken to be identical to the mass of the projectile 
impacting the aluminum target (3*56 X 10"^ slug; 51-95 mg). The dynamic coeffi- 
cient of viscosity of the steel plate is 300 [1.4364 \ and the dynamic 

ft 2 \ cm 2 / 

yield strength 200,000 psi (137*9 -^ 77 V Note that the weight per unit area of 

\ cm 2 / 
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Figure 3.- Radius of perforation as a function of initial projectile velocity for aluminum 
plate, h = 0.25 inch (0.635 cm); Mp = 3*56 X 10" ^ slug (51.95 mg); p = 100 

(o.Vf88 p = 5.2 (2.68 -§3} k = 100,000 psi (68.95 a = 0.*688 inch 

\ cm 2 / ft 5 \ cm 0/ \ cm \ d J 

(0.119 cm) . 


the steel plate is equivalent to the weight per unit area of the aluminum plate. 
Thus, the momentum exchange between the impacting projectile and both the steel 
and aluminum plates is identical. 

As can be seen from figure 6, including material yield strength in the 
solution still has little effect on the perforation radius over the range of 
velocities shown. (The maximum difference between the case K ^ 0 and the 
case K = 0 is of the order of 5 percent.) Although steel has a higher yield 
strength than aluminum, the critical strain rate, which depends on viscosity as 
well as yield strength, is smaller. Hence, the perforation radius in steel is 
larger than in aluminum (for the same initial velocities) and separation occurs 
after a greater passage of time. At the lower velocity range the prevention of 
perforation is again noted as for the aluminum target since for a specific 
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initial projectile velocity of 526 fps (l 60 m/s) the results of the analysis 
indicated no perforation. 



q 0 , 4 • '• r-V r; - 


Figure 6 .- Radius of perforation as a function of initial projectile velocity for steel 
plate, h = 0.088 inch (0.2235 cm); Mp = 3.56 X 10 " 6 slug (51.95 mg) ; p = 300 

1.436 f. = 14.768 — ( 7.61 -S-l k = 200,000 psi 

cm ft* \ cm J) 

a = 0.04688 inch ( 0.119 cm). 

The effect of including yield strength in the analysis when both solutions 
(K = 0 and K ^ 0) yield perforations has been shown to be insignificant. How- 
ever, in the lower velocity ranges the simplified solution for K = 0 may indi- 
cate a perforation whereas the present solution for K f 0 does not. This 
result suggests that the present solution be applied in the case of the determi- 
nation of minimum plate thicknesses necessary to prevent perforation. 


157.90 


kN 


Perforation Radius as a Function of Projectile Mass 
and/or Plate Thickness 

In figure 7 is shown the effect of varying the thickness of the aluminum 
plate or the projectile mass on the nondimensional perforation radius. The 
specific initial velocity of the aluminum projectile was taken to be 20,000 fps 
(6,096 m/s) and the plate thickness was varied from 0.25 inch (0.655 cm) to 


15 



2.50 inch (6.55 cm), all other parameters being held constant. The increase in 
plate thickness decreases the perforation radius and time required for separa- 
tion. Including yield strength in the analysis has little effect on the calcu- 
lated perforation radii (maximum variation in calculated radii is of the order of 
k percent). Figure 7 also represents the effect of decreasing the mass of the 
projectile on the nondime ns ional perforation radius. If the mass of the aluminum 
projectile is decreased from 5*56 X 10 “6 slug (51*95 rag) to 0.556 X 10 slug 
(5.195 mg) (the radius remaining constant) instead of the plate thickness being 
increased from 0*25 inch (0.655 cm) to 2*50 inch (6.55 cm)^ the calculated re- 
sults would have identical plots. This result is due to the fact that the nondi- 

mensional parameter H (ratio of twice the mass of the plug directly beneath the 


projectile to the sum of the mass of the projectile and the mass of the self-same 
plug 2 \ changes in the same proportion with a decrease in Mp or an 

Tta^ph / 



Figure 7.- Effect of variation in aluminum plate thickness 
and/or projectile mass on perforation radius. 

H = 100 slug - (0.4788 6o = 20,000 fpe (6,096 -2-'); 

ft- sec V era 2 / V sec / 

p = 5.2 (2.68 -£-\; k = 100,000 psi (68.95 

«5 V ™ 3 / V cm 7 

a = 0.04688 inch (0.119 cm). 
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increase in h. The perforation radius and separation time are thus seen to 
decrease with decreasing projectile mass. It is evident from figure 7 that the 
effects of yield strength need not he considered in the calculation of perfora- 
tion radii. 


Perforation Radius as a Function of Coefficient 
of Dynamic Viscosity 

The variation in nondimensional perforation radii due to changes in the 
dynamic coefficient of viscosity is shown in figure 8 for the 3-56 X 10“^ slug 
(51*95 m g) projectile impacting on the aluminum plate. A specific velocity of 
3,940 fps (l,201 m/s) was chosen and the dynamic coefficient of viscosity was 



Figure 8.- Effect of variation in coefficient of dynamic viscosity on perforation radius, 
h = 0.25 inch (0.635 cm); Mp = 3.56 x 10"^ slug (51*95 rog); g Q = 3 > 9^0 fp E ^1,201 

p = 5.2 ( 2.68 -£-\; k = 100,000 psi 68.95 — i a = 0.01+688 inch (0.119 cm), 

ft-' \ cnK/ \ cmy 
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varied from 0 to 200 ^~ s . ^ c [0.958 ]; all other parameters remained constant 

ft 2 \ cm 2 / 

and are given in figure 8. Although calculations were made for larger values of 
[i, it was found that the dimensionless separation times were continuing to 
increase above t = l/2; and thus these results were not considered valid 
because of the M short time" nature of the solution. (See appendix.) 

The variation in the coefficient of viscosity is shown to produce a very 
large variation in perforation radii and indicates a definite need for better 
evaluation of this coefficient. The inclusion of the yield strength in the 
analysis exhibits results for perforation radii as much as 10 percent lower 
than the solution neglecting yield strength. However, this variation can be 
considered to be negligible in comparison with the variations of Fp because of 
the lack of accurate knowledge of the coefficient of viscosity (in the present 
state of the art) at these high impact velocities. 


EFFECTS OF TARGET -MATERIAL YIELD STRENGTH ON 
VELOCITIES, DISPLACEMENTS, AND STRESSES 


In the previous sections it has been shown that the results for perforation 
radii obtained with the inclusion of yield strength in the analysis differed 
little from the results obtained when the yield strength was neglected. The 
calculations of the strain and strain rates are little affected by the inclu- 
sion of yield-strength terms and result in similar curves when plotted as a 
function of plate radii. The velocities, displacements, and stresses of the 
solution containing yield strength, however, are at times markedly different 
from the solution in which yield strength is neglected. 


Velocities and Displacements 

In figure 9 are shown the velocities and displacements plotted as a func- 
tion of the plate radius for various times after impact. The pertinent param- 
eters of the problem are those of the aluminum plate discussed earlier and are 
given in figure 9* (The strain and strain-rate distributions are given in 
figures 2 and 3 for K ^ 0.) The specific initial velocity of the projectile 
is taken to be 3^ 9^0 fps ( 1,201 m/s). In figures 9(a) and 9( c ) the results are 
shown for the solution neglecting yield strength whereas figures 9(b) and 9(d) 
show the results obtained by including yield strength. It can be seen that the 
two solutions begin to differ immediately and increase in variance with 
increasing time. The solution neglecting yield strength possesses no rigid 
region since dv/dr never reaches zero (except at infinity), and the material 
remains visco-plastic everywhere. Note also, since cSv/dr is never zero, 
r rz is always greater than the yield stress (which in this case is zero) and 
equation (2b) for r rz is valid everywhere. The solution including yield 
strength, however, develops rigid and visco-plastic regions which are time 
dependent. A Bingham solid becomes rigid when dv/dr vanishes (see ref. 6, 
p. 138 ) and dv/dr vanishes at finite distances from the point of impact when 
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the yield strength is included in the analysis. The second boundary condition 
(eq. (8)) states that as r -*oo, V -> 0. Hence, when dv/dr vanishes at a 
point the velocity nrust remain constant, and for the condition at infinity to 
be satisfied this constant must be zero; therefore, the rigid region possesses 
zero velocity. Consequently, it is concluded that when either the velocity or 
dV/dr vanishes, the material is considered rigid. 

Comparison of the velocities with and without yield strength (figs. 9(a) 
and 9(b)) indicates a much greater velocity decrease with time when the yield 
strength is included. Similarly, the displacements (figs. 9(c) and 9(d)) cal- 
culated from an analysis including yield strength are much less than those 
obtained when yield strength is neglected. 



(b) Velocity as a function of radius for K = 0*551* 
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Figure 9*“ Continued. 



(c) Displacement as a function of radius for K 



(d) Displacement as a function of radius for K = 


Figure 9-- Concluded. 




(a) Velocity as a function of radius for K = 0. 

Figure 10.- Velocity and displacement as a func- 
tion of radius for K = 0 and K = 2.226. 
h = 2.5 inch (6.55 cm); M x = 3-56 x 10~6 slug 

(51.95 mg); P = 100 slug ( 0.4788 

ft-sec \ cn. 2 / 

g 0 = 3 , 9^0 fps (l ,201 |); p = 5.2 

v > ft - 1 

(2.68 — S_\; k = 100,000 psl ( 68.95 —V 
V cm^/ \ cm^/ 

a = 0 . 01.688 Inch ( 0.119 cm). 
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(b) Velocity as a function of radius for K = 2.2; 
Figure 10.- Continued. 



If the projectile mass is decreased to one-tenth its value (or the plate 
thickness is increased tenfold) and the perforation radius is evaluated, it is 
found that the projectile fails to perforate the plate. The results of such a 
decrease in projectile mass (or increase in plate thickness) are shown in fig- 
ure 10 for the aluminum plate. Again, figures 10(a) and 10(c) show results 
obtained by neglecting yield strength whereas figures 10(b) and 10(d) show 
results obtained by including yield strength. 

Note that the solution including yield strength possesses both a visco- 
plastic region and a rigid region, whereas for K = 0, the rigid region has 
retreated to infinity. Also note that when the yield strength is neglected, 
the velocity at ¥ = 1 approaches a limit as time increases but does not drop 
to zero. However, when the yield strength Is included (see fig. 10(b)), the 
velocity drops to zero for all r and the material is considered rigid every- 
where and successfully stops the projectile. The displacements at this point 
(fig. 10(d)) indicate a definite bulge. The effects of yield strength in 
general again indicate much greater velocity decreases and much smaller dis- 
placements than when the yield strength is neglected. 



(c) Displacement as a function of radius for K = 0. 
Figure 10.- Continued. 
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Stresses 


When the radial distribution of shear stresses present in the plate after 
impact are computed from the analysis with and without yield strength, it is 
obvious from equation (2b) that they will differ. In fact, the variance at the 
time and radius of perforation is approximately 2 to 1 . In figure 11 is shown 
a distribution of shear stress computed at the dimensionless time of t = 0.25 
(calculated separation times are actually t = 0.24 for K = 0 and t = 0.23 
for K / 0) for a l/4 inch (0.635 cm) aluminum plate. The initial projectile 
velocity was taken to be 3,9^0 fps ( 1,201 m/s) and the radius of perforation at 
separation is indicated by the dashed lines in figure 11 . As can be seen from 
figure 11, the shearing stresses occurring at the radius of perforation computed 
by use of the two different solutions differ by a ratio of about 2 to 1. 


CONCLUDING REMARKS 


It has been shown in the present paper that the inclusion of target - 
material yield strength in a one -dimensional analysis of hypervelocity impact 
perforations has produced little effect on the resulting perforation radii. 

The variation in calculated perforation radii as compared with the solution 
neglecting target yield strength amounts to only 5 percent. This percent dif- 
ference in calculated perforation radii remains at the negligible 5~P erc ent 
level even with variations in the pertinent parameters such as plate (or target) 
thickness, projectile mass, and initial velocity. With variations in the 
dynamic coefficient of viscosity, the difference does increase to perhaps 
10 percent but the variations in perforation radius due to differences in the 
assumed value of the viscosity coefficient alone are much greater than the 
10 percent present when yield strength is neglected. In fact, the determina- 
tion of accurate values for the dynamic coefficient of viscosity is much more 
critical in the calculation of perforation radii than is the inclusion of the 
target yield strength. 

In the determination of the minimum thickness of plate necessary for pre- 
vention of perforation, however, the effect of yield strength may be important 
depending on the initial velocity and mass of the projectile, and the present 
analysis containing yield strength should be applied. Furthermore, parametric 
studies in which not only the initial projectile velocity is varied but also the 
thickness of the plate and the mass of the projectile could provide pertinent 
data as to minimum thicknesses of target material needed to prevent perforations. 

The velocities, displacements, and stresses can be markedly different in 
the resulting calculations for hypervelocity impact depending upon whether 
target yield strength has been included in the analysis or not. The effects 
of yield strength in general are to produce much greater velocity decreases 
and much smaller displacements than those obtained when the yield strength 
is neglected. In fact, for the particular case where the projectiles do not 
perforate the target plate, the analysis including yield strength indicates 
that the velocity drops to zero and the plate is deformed, whereas if the yield 
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strength is neglected, the velocity only approaches a limiting value, which is 
not zero. Thus, for this particular case, the velocities, displacements, and 
stresses calculated hy using the simplified solution may be grossly in error. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va. , September 30, 


1964 . 
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APPENDIX 


GOVERNING EQUATION OF SHEAR PERFORATION 


An analytical solution is presented for the governing linear differential 
equation and its associated houndary conditions. In a manner similar to that 
employed by Chou in reference 3, the governing partial differential equation is 
reduced to a total differential equation by use of Laplace transform techniques. 
To obtain a particular solution to the resulting total differential equation, a 
further change in variables becomes necessary. After satisfying the associated 
boundary conditions, a "short time" solution is determined by employing asymp- 
totic approximations in the transformed state. 

The governing linear differential equation of motion in the axial direc- 
tion, for simple shear perforation, as derived in 
written as 

c^V 1 dV 1 dV = 1 

£ r 2 r Sr “ v St r 

Transforming equation (Al) with respect to t by 
techniques results in 

+ 1 a£ . l[ s v - V(0,r)] 
aJ2. r dr v <- u 


where s is the transform parameter. By use of the initial condition 
(eq. (4)) V(0,r) =0 at t = 0, r > a, equation (A2) is reduced to 

d 2 V , 1 dV sV _ 1 Jl (A3) 

dr 2 r dr V r ps 


the text (see eq. (3)) can be 


k 


(Al) 


use of Laplace transform 


1 _k_ 
r (is 


(A2) 


Let 


r = 



(A4) 


Hence 



APPENDIX 


a 2 ? __ _afv s 
dr^ dz2 v 

and by substitution of these relations into equation (Aj) there results 


(A 6) 


+ V = (A7) 

dz 2 z dz £ |ISZ 


The homogeneous solution for V(z) can he written directly as 


V H (z) = AJq ( z) + BY 0 (z) 


(A8) 


For the particular solution to equation (A7) the Struve function of zero order 
is satisfactory. A known relation for the Struve function of order p is 
(see ref. 7; P* 211) 


d 2 S p (z) x dSp(z) 

— — — + =■ -4- — + l 


dz" 


z dz 


- £?)%<’ 


) - 


Ml 


iRp - 0.5): 


(A9) 


Let p = 0; then 


d 2 S 0 (z) 2 . dS Q (z) 


dz" 


+ — + Sn(z) = 

z dz uv z It 


(A10) 


Hence, the particular solution to equation (A7) for V(z) can be written as 


V-n(z) = 


^o( z ) 

[IS 


ft 


(All) 


and the general solution becomes 


V(z) = V H (z) + V p (z) = AJ 0 (z) + BY 0 (z) 2p_S 0 (z) 




or resubstituting for z from equation (a4) 


(A12) 


V(r,s) . AJ 0 (ir /|) + EY 0 (ir /l) - ^M 0 (ir/I) 


(A13) 
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In order to apply the (transformed) boundary condition given by equa- 
tion ( 8 ), namely, V = 0 as r -» <», one must consider the asymptotic behavior 
of the solution as given by equation (A13) for |r| » 1. The asymptotic 
approximation for the Struve function of order zero as found in equation ( 136 a) 
of reference 8 is 


S o( lr 0 ) - *o( lr 0 ) 


* lr it 


(Al4) 


and since a known identity between the Bessel functions is (see eq. (.110) of 
ref. 8 ) 


then 


K °( r J!) - ffoC 1 - /I) + “4 r if) 


S o( lr if) ” i:t o( r Jf) - fro(r if) 


(A15) 


*ir |f 


(A16) 


Hence, for | r | » 1, the general solution can be written as 
V(r,s) « AI 0 (r 


if) + B [ lI o( r if) - ¥ °( r if) 


rtk r.L |1 

3 \ . ik _ir 7 r f 

in i i - 

rr U V 

'j + 1 7 °V v 

2||xs ^ 

l^ s 


s 2 

r^i 

v 


(A17) 


Application of the transformed boundary condition V = 0 as r -»<» (see 
eq. ( 8 )) yields the relation 


A = -Bi - 


otk 


2 0 s 


(Al 8 ) 


since rvQ — > w anu. iq — > w as ± — * ^ . use wi une xviiewii ^ ^ ^ 

Ko(r |f) (eq- (A15)), s o( ir j|) (eq- (Al4)), I 0 ( r J|), and Y 0 ( ir jf) 


since Kq -> 0 and Iq -> 00 as r -»». Use of the known identities between 

«o(r (eq- (A15)), 

reduces the general solution to 

V(r,s) = -B% 0 ( r ||) " 


irk 


2 if* s 


(r Jf) - Jf) C«9) 


2 if 13 
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Application of the (transformed) boundary condition (see eq. ( 7 ); note also 

aq- (5)) 

sV - fu— + t— = V 0 (r = a, t > 0) (A20) 

dr s u 

permits the evaluation of the quantity B. Substitution of equation (A19) into 
equation (A20) and use of equation (l40a) of reference 8 yield the following 
expression for B: 


B 



(A21) 


Hence, the general solution can he written as 


V(r,s) 


f * ^°(* $ * s °(‘ a ^ ^ ® + E ® 

1 - Sifia 
Tt ± \ V v J\ 

"3 

[4(f) 

sK 0^ a \|f^ + r > /I 

) 





jtik 



(A22) 


A "short time" solution is now determined from the general solution by 
assuming a ^ >1. By reference again to equation ( 156 a) of reference 8 , it 

can be shown that 


S l(- |l)” Y l( ia #) 




(A23) 


A known identity between the modified Bessel function of the second kind of 
order one Kq and the Hankel function of the first kind of order one H-^ is 
(see eqs. ( 115 ) and ( 70 ) of ref. 8) 


K i( a if) * - ( la if) - -f[ J i ( Ia if) + 1Y i ( la if) 


(A24) 
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Hence 

J 1 ( ia #) ' '“I ( ia It ) - ^1 ( a /!) < A2 5) 

and since 

then, from equation (A25) 

Ix(a If) - l/l) + /I) (A27 > 

Thus, 

| - Sl (la H) - Il(a /I) ♦ (a |f) («8) 

2 

Substitution of the asymptotic expressions for S q and — - S]_ into the 

general solution (eq. (A22)) and making use of the identity for Yq (eq. (A15)) 
yield 


V(r,s) 



(A29) 


Note that for increased accuracy in equation (A29), it is only necessary 
to increase the number of terms taken in the asymptotic expansions of Sq and 
Si; of course, the semi convergent nature of the asymptotic expansions must be 
considered. These additional terms are polynomials and would not present any 
new difficulties. To aid in determining the inverse transform of equa- 
tion (A29), the following asymptotic expansions of the modified Bessel func- 
tions Kq, Ki are employed (see eq. (ll4) of ref. 8) 


Kq(z) 

- |jLe~ z | 

V2z 

fi - r + 

9 

V 8z 

128z 2 

K x (z) 

-z 

V2z^ 

(i + ± - 

15 

V 8z 

128z 2 



•) 


(A50) 
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Substitution of the asymptotic expansions given in equation (A30) into equa- 
tion (A29) yields after algebraic manipulation: 


V(r,s) 



V 2 



+ 



- 2ar - 7r 2 
128a 2 r^ 


+ 


a - 3r 
8ra 



1 

s_ 

V 


+ 


k 

ps^r 


(A31) 


and by application of standard transform tables (ref. P* formulas ( 11 ) ) 

V becomes 


V(r,t) 


= V, 



1/2 


erfc 




erfc 

b v J a 


r - a 
2<fvt 


9a 2 - 2ar - Jr^ 
< 128 r 2 


Ua a - 

^ v 8r ^ 


H 2 a^ 


^i 2 erfc^ 


a^ 


2 ft 


'i. - rV*7 /2 




/ 9a^ - 2ar - 7 r ‘ 

+ I — + 


128 r 2 


4ti 2 erfc + fe-i® - ^jE(Ut) 5 / 2 i 5 erfc 

L 2 fi \ 8r v / a 2^vt 

.jxa a_-_5r 2 4 L-L± . . . 

8r v 2 /^ 2 Jft 


k|i /aW 2 


l6t 2 Perfc r 


2 Jvt 


FT + (^ - ^ (4t) 


5/2,5, 




+ ( 9a 2 - 2ar - 7r 2 + ^a a_-_5r + ^2 ^Y^tSArfc 
\ 128 r 2 v 8r v 2 /a 2 


i^erfc 


r - a 


2\jvt 


r - a 
2 \fvt 


kt 

pr 


((r - a) > o) (A32) 

where the symbol i n denotes the nth integral of the complementary error 
function. 
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Equation (A32) can now be nondimensionalized as follows: 


Let 


— _ r 
a 


T* _ Vt 


tt *ua _ 2ata^hp 

H = ^ = “IT” 




(A35) 


K = 


V 

ka 

V 


J 


Hence, the nondimen sional velocity can be written as 


J-.U 


2)ffi 


erf c + tK(l - H)ti 2 erfc r ~ - + l6KHt2i^erfc 


+ Ri 2\/ti erfc 


2)jf 

r ' ■■ I + 8 k(i - H)t5/2i5 


ft 


erfc - ~_ X + 32KHt 5 / 2 i 5 erfc 


r - 1 


>.w 


ift. 


+ R 2 [4ti 2 erfc r ~_ X + l6K(l - H)t 2 i 4 erfc r ~ — + 64KHt3j6 er fc r ' 1 


L 


2ft 


2ft 


2<[t 

((r -1) 2 0) 


Kt 

F 


(A3*0 


where 


R i ■ 5 - i - H 


r 2 = 


128?^ 64 r 128 


+ H 


(M) 


+ 


By comparison of equations (A3l) and (A 32 ) } it can be shown that s is related 

to l/4t and that a\[sjv > 1 implies t < l/4. Nevertheless ^ calculations 
indicate that the three-term asymptotic series presented in equation (A30) are 

only approximately 10 to 15 percent in error for t = l/2 ^and a\fs/ v = l/\/2^. 
Therefore, results are presented in this paper for values of t as high as l/2. 
Results for t > 1/2 should be viewed with increasing skepticism. Differentia- 
tion of equation (A3^0 with respect to r yields the nondimensional shear 
strain rate 
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dv &_ _ _1_ 
v ° " 


(r-l)2 



+ 


erf c 


r “J- + 2 ft |R4 - K(l - h[| i erfc 
2ft 



+ 4t[R 5 + K(l 


H)R^ji 2 erfc — 


-r- + 8(t)5/2[ K( i _ h)R 4 - KHji5 erfc 


r - 1 






+ l^fKHR, + K(1 - H)R.n i^erfc £ ■ '■ 1 + 32(t) 5 / 2 KHR4i 5 erfc £-~J- 
L 2^t 


2\ft 


i^erfc r ~ ^ . . A + — 

5 2jf ( ^ 


((r -l)io) (A55) 


where 


r 5=-J;-5 + h 


r 4 = 


_ 15 


5 + J- + 


128T 2 64r 128 


Sl*’)- 




Rr = 


*+5 + 3 + _7_ + _h/5 _ _3_ 

256r 5 128? 2 2 ^6r 16 \r T 2 


h£ 

2r 


By integration of equation (A35) with respect to t, the nondimensional shear 
strain is obtained 


dw y 
dr V 0 a 


= iilA-L 
\|? ^ 2\fF 


■i erfc 


r + R 5 i 2 erfc r "j 1 + 2||t (r 4 - K(l - hH i 5 erfc 
2\|t~ 2{t 2\ft 


+ 


4tjR 5 + K(l - H)R^] Aerf c 


£-=-i + 8(t)5/ 2 jK(l - H)R 4 - KHj i^erfc Z-zJL 

2$t 2jf 


+ l6t 2 [KHR 3 + K(l - H)RrJ i 8 erf c r - + 32(t) 5 / 2 KHR 4 i 7 erfc r ~- l 

2 ft 2\jf 


+ 64t?KHR 5 i 8 erfc 

2W 


k£ 

2F 2 


((F - l) ^ 0 ) (A36) 
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A second integration with respect to r would determine the corresponding dis- 
placement. However, a simpler approach is to integrate the nondimen sional 
velocity with respect to t to obtain the nondimensional displacement 


r v = hi /j2 er f C r — + 4K(l - H)ti 4 erfc + l6KHt2i 6 erf c £ ~~ 

a 2v o 2\ft" 

+ RXjtlWc + 8K(1 - H)(t) 5 / 2 i 5 erfc + 32 KH(t) 5 / 2 i 7 erfc 

L L 2^ 2^ 2^- 

+ R 2 [tti 4 erfc + l6K(l - Hjt^erfc £-=^ + 64KHpi 8 erfc ■ • 

L 2 ft 2jt 2n J 

((r - 1) £ 0) (A37) 



The symbol i n erfc x represents the repeated nth integrals of the complementary 
error function erfc £ 


i n erfc x = f i n ~^erfc £ (A 58 ) 


in which 

i°erfc x = erfc x 


i _1 erfc x = -^e x 

and x(|) is the argument of the complementary error function. A recurrence 
relation useful in determining the repeated integrals of the complementary error 
function is (see refs. 8 and 9) 


i n erf c x = i n ~ 2 erfc x - 2xi n ^rfc x 

2n 


(A39) 


A table of the error function and its derivatives and integrals for values of 
the argument between 0 and 3-0 is presented in reference 9 from values given in 
reference 10. A similar table has been included in this paper, with slight 
extensions, for convenience in computing these expressions. (See table II.) A 
table of the first 11 repeated integrals of the error function with values of 
the argument between 0 and 2.80 is also available. (See ref. 11*) 
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TABLE II.- FIRST EIGHT REPEATED INTEGRA IB OF THE COMPLEMENTARY ERROR FUNCTION 


X 

erf ' x 

erfc x 

i erfc x 

i 2 erfc x 

i^erfc x 

i4erfc x 

i5erfc x 

i^erfc x 

i^erfc x 

i^erfc x 

0.00 

1.12837910 

1.00000000 

0.56418955 

0.25000000 

0. 09403159 

0.03125000 

0.00940315 

0. 00260416 

0.00067165 

0.00016276 

.05 

1.12556170 

.94362802 

.51559945 

. 22301702 

.08221629 

.02684942 

.00795313 

. 00217117 

.00055257 

. 00013224 

.10 

1.11715160 

.88753708 

.46982210 

.19839317 

.07169057 

.02300688 

. 00670892 

.00180542 

.00045341 

. 00010717 

.15 

1.10327410 

.83200403 

.42683645 

. 17598827 

.06233999 

.01966078 

. 00564417 

. 00149729 

.00037107 

. 00008662 

.20 

1.08413470 

.77729741 

. 38660787 

.15566357 

. 05405707 

. 01675509 

, 00473550 

.00123840 

.000*0286 

. 00006982 

.25 

1.06001410 

.72367361 

. 34908865 

. 13728232 

. 04674124 

.01423896 

.00396217 

.00102148 

.00024653 

.00005613 

■ 50 

1.03126090 

.6 713752+ 

. 31421848 

. 12071054 

.04029869 

. 01206641 

.00330588 

.00084024 

. 00020012 

.0000450I 

.35 

.99828371 

.62061795 

. 28192557 

.10581751 

.03464221 

. 01019599 

.00275050 

. 00068922 

. 00016200 

.00003598 

.40 

.96154130 

.57160765 

. 25212759 

. 09247640 

. 02969107 

. 00859044 

. 00228187 

.00056374 

.00015077 

.00002869 

■ ^5 

.92153201 

. 52451828 

. 22473278 

. 08056470 

.02537075 

.00721637 

. OOI8876O 

. 00045979 

. 00010527 

.00002281 

.50 

.87878258 

. 47950012 

.19964123 

.06996473 

. 02161275 

.00604399 

. 00155687 

.00037392 

. 00008449 

.00001808 


.83383665 

.43667663 

. 17674618 

. 06056396 

.01835430 

.00504677 

. 00128028 

.00050320 

.00006762 

.00001430 

.60 

.78724343 

. 3961+391 

•15593557 

.05225536 

.01553815 

. 00420119 

.00104967 

. 00024513 

.00005396 

.00001127 

.65 

.7595^676 

• 35797067 

. 13709245 

.04493762 

. 01311225 

.00348646 

. OOO85798 

. 00019758 

. 00004293 

.00000886 

.70 

. 69127486 

. 32219881 

. 12009827 

.03851530 

. 01102947 

. 00288425 

. 00069915 

. 00015878 

.00003406 

. 00000694 

■75 

.64293107 

.28884366 

. 10483279 

. 03289861 

. 00924747 

. 00237842 

.OOO56798 

. 00012720 

. 00002694 

. 00000542 

.80 

.59498579 

. 25789904 

.09117366 

.02800529 

.00772753 

.00193515 

. 00045992 

. 00010160 

.00002125 

. 00000422 

.85 

. 54786972 

.22933194 

.07900272 

.02375682 

.00643601 

. 00160194 

.00037127 

. 00008089 

. 00001669 

.00000328 

.90 

. 50196857 

.20309179 

.06820167 

. 02008219 

.00534228 

.00130826 

. OOO29874 

. 00006421 

.00001308 

. 00000254 

■ 95 

. 45761925 

. 17910919 

.05865589 

.01691575 

.00441932 

. 00106487 

.OOO2396O 

. 00005080 

.00001022 

. 00000196 

1. 00 

.41510750 

. 15729921 

.05025454 

.01419753 

.00364324 

.OOO86388 

.00019154 

. 00004006 

.00000795 

.00000150 

1.10 

• 33647960 

.11979493 

.03646538 

.00989277 

.00245021 

.00056278 

. 00012120 

. 00002467 

.00000477 

.00000088 

1.20 

.26734435 

. 08968602 

.02604895 

.00679213 

.00162463 

.00036162 

. 00007567 

. 00001500 

.00000283 

. 00000051 

1.50 

.20820799 

. 06599206 

.01831432 

.00459370 

. OOIO6178 

.00022913 

. 00004660 

. 00000899 

.0000016;} 

. 00000029 

1.40 

. 15894171 

. 04771488 

. 01267002 

.00305970 

. 00068380 

.00014312 

. 00002830 

.00000532 

.00000095 

.00000016 

1.30 

.11893029 

.03389485 

.00862287 

. 00200656 

.00043386 

. 00008812 

. OOOOI695 

. 00000310 

. 00000054 

. 00000009 

1.60 

. 08722906 

.02365162 

.00577193 

.00129535 

.00027113 

.00005346 

. 00001000 

. 00000178 

.00000030 

.00000005 

1.70 

.06271104 

.01620954 

.00379930 

. OOO82297 

. 00016686 

[ .00003195 

. OOOOO582 

.00000101 

. 00000016 

.00000002 

1.80 

. 04419172 

. 01090950 

. 00245876 

. 00051449 

i .00010109 

. 00001881 

. 00000333 

.00000056 

. 00000009 

.00000001 

1.90 

.03052474 

.00720957 

. 00156418 

.00031641 

! . 00006030 

. 00001090 

. 00000188 

. 00000031 

. 00000004 

.00000000 

2.00 

.02066699 

. 00467774 

.00097801 

.00019142 

i .00003538 

.00000623 

.00000104 

. 00000017 

. 00000002 

. 00000000 

2.10 

.01371563 

. 00297947 

.00060093 

.00011388 

: .00002043 

.00000350 

. 00000057 

. 00000009 

. 00000001 

.00000000 

2.20 

. 00892216 

. OOI86285 

.00036281 

.00006662 

: .00001161 

. 00000194 

.00000030 

. 00000004 

. 00000000 

.00000000 

2.30 

. 00568902 

.00114318 

. 00021519 

.00003831 

j , 00000648 

. 00000105 

.00000016 

. 00000002 

.00000000 

. 00000000 

2.40 

.00355565 

.00068851 

. 00012540 

.00002164 

.00000358 

. OOOOOO55 

.00000009 

. 00000000 

.00000000 

-.00000000 

2.30 

.00217828 

.00040695 

.00007176 

. 00001203 

.00000193 

. 00000029 

. 00000004 

. 00000000 

. 00000000 

-.00000000 

2.60 

.00130805 

.00023603 

. 00004034 

. 00000655 

. 00000104 

.00000014 

. 00000003 

-.00000000 

. 00000000 

-.00000000 

2.70 

. 00076992 

. 000134-33 

.00002226 

.00000351 

. 00000054 

.00000007 

. 00000001 

- . 00000000 

. 00000000 

- . 00000000 

2.80 

.00044421 

. 00007501 

. 00001207 

. 00000184 

. 00000029 

.00000002 

. 00000001 

-.00000000 

. 00000000 

-.00000000 

2.90 

.00025121 

. 00004110 

. 00000641 

. OOOOOO97 

.00000012 

.00000002 

-.00000000 

. 00000000 

-.00000000 

.00000000 

3-00 

.00013925 

. 00002209 

. 00000335 

. 00000049 

.00000006 

. 00000000 

. 00000000 

.00000000 

. 00000000 

-.00000000 
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